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Abstract 

The multifractal spectrum of various three-dimensional representations of Packed Swiss 
Cheese cosmologies in open, closed, and flat spaces are measured, and it is determined that 
the curvature of the space does not alter the associated fractal structure. These results are 
compared to observational data and simulated models of large scale galaxy clustering, to 
assess the viability of the PSC as a candidate for such structure formation. It is found that 
the PSC dimension spectra do not match those of observation, and possible solutions to 
this discrepancy are offered, including accounting for potential luminosity biasing effects. 
Various random and uniform sets are also analyzed to provide insight into the meaning of 
the multifractal spectrum as it relates to the observed scaling behaviors. 
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1 Do We Live in a Fractal Universe? 



The notion of a hierarchically-structured world is a recurrent theme in our 
understanding of Nature |T]. According to the Cosmological Principle, the 
Universe must be homogeneous and isotropic. This oft cited-as-fact stip- 
ulation is the basis for the Friedmann-Robertson- Walker (FRW) solutions 
to Einstein's Field Equations, from which the expansion dynamics of the 
Universe are derived [2|. However, the breakdown of homogeneity (at least 
on smaller distance scales in the Universe) is quickly becoming an accepted 
ideal in cosmological circles. That is, while homogeneity requires matter to 
scale uniformly in space (i.e. Dp = 3, equal probability scaling behavior 
in all spatial directions), actual measurements of the distributions suggest 
otherwise. 

The existence of fract ally-clustered matter was further emphasized by 
Peebles (see [3] ) , whose various two-point correlation analyses of three-dimensional 
catalogs yielded the exponent 7 ~ 1.7, the co-dimension of which was taken 
to be the fractal dimension Dp = 3 — 7 ~ 1.3. Some successive works 
confirmed this value, showing a study of the CfA redshift survey to agree 
statistically |U|3]. However, these same studies cite model simulations with 
higher dimensionality, in particular D ~ 2 1 . 

Over the course of a decade, increasing evidence has been put forth to 
suggest that the dimensionality of a wide range of redshift survey catalogs 
yield such clustering dimensions. In the comprehensive publication j^j, the 
results of fractal dimension analysis of numerous galactic catalogs are re- 
viewed, with the general consensus that each data set reveals a unanimous 
Dp ~ 2 (see Table EJ. 

Since each catalog is limited in size and spatial extent (i.e. volume), the 
associated dimensions can only be statistically viable up to some effective ra- 
dius R e s (the radius of the largest sphere one can inscribe in the associated 
sample, without surpassing the catalog boundaries). For most galaxy cata- 
logs considered in the literature, these range up to 50 h~ l Mpc (recall that 
h^ 1 is the scaling factor in the Hubble Constant, H = 100 h km s _1 Mpc _1 
jU]). Pietronero et al. further extend the analysis to include clusters and 
superclusters, concluding that the two sets of data represent the same self- 

1 This discrepancy ultimately has its roots in the method of calculation. The dimension 
D = Dp = 2 is the standard fractal dimension obtained by usual means, while The value 
D r*> 1.3, calculated from a two-point correlation function, is actually the q = 2 multifractal 
dimension D-i- See Section 0] 
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similar structure at differing scales. While clusters are more distant, they 
are more luminous, and by counting the clusters as single objects, a coarser 
estimate of scaling is obtained over larger distances. This provides evidence 
for fractal behavior up to distances 1000 h^ 1 Mpc, with no suggestion that 
homogeneity ensues. 

Of course, a potential limitation of such estimates is the lack of a definite 
third dimension. While angular spans may be accurately measured (bar- 
ring external interference from sources such as gravitational lensing), the 
(co-moving) distance to the objects in question must rely on estimates from 
luminosity-distance relations or redshift measurements. These effects are in- 
troduced primarily by the recessional velocities of the objects, and in general 
may be obtained via application of the Hubble redshift law, whose simplest 
form for low redshifts is linear, 



and for general redshifts (Mattig's 1958 relation; see e.g. jH]) 



in an expanding Universe with deceleration parameter q . For q = 1/2 (i.e. 
flat Universe), this can be shown to reduce to |Hj 



Typical redshifts for the majority of nearby catalogs are z < 0.05, but deeper 
surveys such as ESP or LCRS can contain redshifts of the order z ~ 0.2 [7JIH]. 
Note that there is no correction for peculiar velocities in the associated figures 
of Table EE]©. 

The authors of [Hj conclude that use of the Euclidean Hubble relation (0) 
in catalogs where the redshifts are high instead of © or (j3J) does not affect 
the estimated dimensionality, and thus the result Dp ~ 2 is stable to such 
variations. This claim is further elucidated in reference 0, who note that 
use of the Euclidean Hubble distance formula may be exported to distance 
greater than 600 Mpc with no consideration of relativistic curvature effects 2 . 

In fact, this provides alternative support for the notion of "deformation independence" 
of the fractal dimension with respect to the embedding manifold, discussed in 



c zq + (q - l)(y/2zq + 1-1) 



(2) 



r = 



Ho q 2 o(l + z) 
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Before proceeding, however, it should be emphasized that the aforemen- 
tioned Dp ~ 2 fractal structure of the local Universe (with no transition to 
Dp = 3) is a hotly debated subject, and is by no means to be interpreted as 
"fact" . Rather, it is a movement largely spearheaded by the authors of [H] , 
whom collectively have published over 50 related articles in the past 5 years. 
In fact, several earlier works [3] suggest that while the local structure may be 
fractal, the transition to homogeneity is clearly marked at about 5 h~ x Mpc. 
The general consensus tends to reflect this finding [i.e. that at least the 
local clustering structure can be described as a fractal), although the transi- 
tion distance varies from catalog to catalog and analysis. More recently, the 
authors of reference ^JJ submit that there is no statistically or physically 
viable way to obtain a fractal scaling behavior for the ESP redshift survey 
and Abell clusters, instead citing a homogeneous D ~ 3 scaling, attributing 
the potential Dp ~ 2 to various redshift-related distance biases on the part 
of Pietronero et ai, including the substitution of the Euclidean Hubble law 
in lieu of the FRW version! A rebuttal to this publication suggests that the 
cross-over in question is explicitly dependent on the cosmological model used 
for the calculations ^2]- Clearly, there is much disagreement in the litera- 
ture, so one must be careful not to take each report at unquestioning face 
value. 

As newer and more comprehensive redshift data becomes available from 
such current surveys as the Sloan Digital Sky Survey (SDSS) [T3] or the 
2dFGRS ^1] (which are targeting well in excess of 100 000 redshifts, and 
in the case of SDSS, 1 000 000 galaxies), the crux of this debate may be 
addressed with more certainty (or, on the other hand, such may serve to 
further complicate the issue!). 

2 The Packed Swiss Cheese Cosmological Model 

The notion of a fractal Universe defies the Cosmological Principle, which 
demands homogeneity and isotropy at all points. Locally, however, it is 
rather evident that the universe is not homogeneous. Any cosmological model 
which is used to represent the observed galaxy distribution must adhere to 
this point, lest its power of predictability be diminished. 

Some of the earlier references to locally inhomogeneous cosmological mod- 
els date back to Einstein and Strauss as well as Schucking in the 1950s 
[To] , and later Rees and Sciama in the 1960s [TS]. Dubbed "Swiss Cheese" 



4 



models, these constitute locally inhomogeneous but globally homogeneous 
spacetimes which everywhere satisfy Einstein's Field Equations. Further 
works studied the effects of multiple hole solutions vis-a-vis gravitational 
lensing effects with and without a Cosmological Constant [THl I17j. and most 
recently the notion of optimally packed, volume-filling hole solutions, dubbed 
Packed Swiss Cheese cosmologies (PSC) [T%1 ITU]. 

A sphere is inscribed in a zero-pressure, expanding FRW Universe of 
spatially-uniform density profile po, and the mass within is condensed to a 
smaller sphere of larger average density p\ > po- The total mass is conserved 
within the shell. Outside the inscribed radius, the space is still purely FRW, 
with line element 



ds 2 = dt 2 -R 2 (t) duo 2 + Sl(oj)dQ 



dt 2 - R 2 (t) 



1 — kr 2 



(4) 



for uj the angular radial coordinate, and the co-moving surface defined by 
Su{oj) is either sin(a;), u, sinh(u;) for positive, flat, or negative curvatures 
spaces respectively (equivalently one can write Sk(oj) = sm(\/kuj)/Vk for 
k = +1, 0, -1) 0HZI. dQ = d6 2 + sin 2 (0) d<p 2 is the standard (spherical) solid 
angle element. Hence, the CP is preserved on the exterior. 

The metric of the vacuum interior of the hole is the Schwarzschild line 
element, 

, 9 / 2Gm Ar 2 \ , 9 dr 2 9 ~ , . 

ds 2 = 1 )dt 2 — — — - r 2 dQ 2 , 5 

V r 3 J 1 - 2Gm/r - Ar 2 /3 v ; 

for general Cosmological Constant A and interior mass m (subject to ap- 
propriate matching conditions at the surface boundary |TB]). If the internal 
mass density has some spatial extent, this can also possess an FRW line el- 
ement (but not necessarily), although in the case of the PSC one considers 
only an interior point mass. Since their average density profile is unchanged, 
any number of holes may evolve independently of each other, since they have 
no gravitational influence on one another (provided the inscribed surfaces do 
not overlap) [T7J. 

The PSC is conceptually similar to the classic "Apollonian Packing" prob- 
lem of efficiently filling a region with tangential circles of varying sizes. In 
the former case, however, it is a three-dimensional packing problem with 
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spheres in spaces of constant curvature. A point is chosen at random in a 
region of FRW space, as before, and a sphere of an arbitrarily large radius 
Ro is inscribed. The mass within this sphere is shrunken down to a smaller 
radius R' , increasing the local density and creating a density discontinuity at 
the boundary. A second point is chosen in the remaining continuous density 
region, and a second sphere is inscribed, subject to the constraint that it be 
tangential to the first. The interior mass is again contracted to a specified 
radius, and the process is repeated indefinitely. As mentioned, the contracted 
radius R' — > 0. 

Any number of packings and configurations can be obtained by varying 
the initial size and placement of the first sphere. Since there are no explicit 
scale constraints on the configurations, a fully recursive PSC can be formed 
by inserting within any sphere of one packing the contents of another. Such a 
model hearkens of self-similarity, at the very least on a statistical level. The 
libraries used for this study (see [Tj3] for details) contain on average between 
30 000 to 90 000 spheres (generally much larger than the size of the redshift 
surveys). Figure Q shows a two-dimensional projection of a typical packing 
library. 

There is no discussion herein of effective scales or cutoff radii, as with the 
survey catalogs, since the packing libraries are scaled in dimensionless units. 
Indeed, since one may recursively pack them at will, the overall physical 
scales may be set arbitrarily. 

It should be noted that this Swiss Cheese model is different than that of 
Ribeiro [2011211122], who also predicts a large-scale fractal structure. Com- 
parisons with the aforementioned results are discussed in section 

2.1 Justification for a Swiss Cheese Cosmology 

The choice of a Swiss Cheese cosmology may seem at first to be without clear 
motivation, but it is in fact one of the most logical and consistent spacetimes 
in which to work. The Swiss Cheese solution is an exact solution of the 
Einstein Field equations at all the length scales involved. This is absolutely 
necessary to begin consideration of fractal behavior since the very essence of 
the fractal approach is validity over a range of length scales. 

The PSC model considered herein is built upon the Swiss Cheese solution, 
and while this solution at first appears to some to be artificial, it covers 
the essential field regions for an object embedded in a background universe. 
Thus it encompasses a near field, a medium field (i.e. the vacuum around the 
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central object) and the far field, which just becomes the FRW background 
universe. The sharpness of the boundaries have little physical impact. Note 
that the sharp boundary of the Earth does not have any drastic effects when 
scene in terms of the gravitational potential (and thus the metric tensor 
either). 

2.2 Cosmological Principle, GR Style 

It is useful to re-formulate the Cosmological Principle in the language of 
General Relativity, insofar as the PSC models are concerned. Recall that the 
curvature of a manifold with metric g a b is defined via second derivatives of 
the metric, and this information is completely contained within the Riemann 
Tensor, 



where [ ] denotes the anti-symmetric sum with respect to index permutation. 
Since it can be shown that C^ ad = 0, the Weyl Tensor is the trace-free portion 
of the Riemann Tensor, and complements the information contained within 
the Ricci Tensor (trace of Riemann). Combined, both objects contain the 
complete curvature information of the Riemann tensor, with Ricci represent- 
ing the local curvature contributions (via Einsetin's equations, i.e. as a func- 
tion of the local pressure p and matter density p via the energy-momentum 
tensor T a b = pg a b + (p + p)u a Ub, with u a the associated four- velocity), and 
Weyl the non-local or external contributions. 

A different interpretation is as follows. For a sphere of radius r inscribed 
in a density field p(r), the Weyl Tensor is a measure of the over /under- density 
of integrated (or enclosed) mass m(r) to "average contained mass" m(r) That 
is, C a bcd oc m(R) — fh(R), where 




(6) 



C a bcd — Rabcd H ~ \9a\dRc\b + gb[ C Rd]a) + 




(n- l)(n-2) 



2 



R 9a[ c gd]b • (7) 



m(R) 
m{R) 




(4n/3)R 3 p(R) , 



(8) 
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and Cabcd is evaluated on some closed surface which encloses the density 
field p (see reference jTH] for a complete derivation). If this quantity is non- 
zero, it signals an imbalance in the mass distribution, and hence provides 
a "gravitational compass". That is, by providing a preferred direction, it 
breaks the isotropy of the space. If, on the other hand, C a ibj vanishes on this 
surface, then the space contains no local tidal forces (no compass), and the 
Cosmological Principle is upheld if the mass distribution is homogeneous. 

Every inscribed sphere in the PSC packings is itself such a surface (see 
Figured}, an d thus it is always possible to find a compound surface within the 
packing on which the Weyl Tensor vanishes. The space outside the inscribed 
boundary is FRW, by design, but that contained within the sphere is not. 
In the PSC model, it is assumed to be a vacuum, and thus has vanishing T ab 
(and thus vanishing Ricci Tensor). 

3 Measuring the Multifractal Spectrum of the 
PSC Models 

Effective scaling dimensions are obtained by the (three-dimensional) box 
counting method. Tables I2HH show the calculated dimension D q , q = 0, 2 (box 
and correlation) and q — > oo, for several packings libraries of flat, positive, 
and negative curvatures respectively. Values of D q for q < are discussed in 
Sectional 

Covering cubes of side d ranged roughly over two orders of magnitude, 
requiring 8 cubes at the largest scales, to about half the population size of the 
packings (i.e. about 1-2 particles per box on average). Below this limit, the 
box scales drop below the average interparticle distance, and the calculation 
becomes skewed by the finiteness of the data set. Since the libraries are 
themselves spherical distributions of points, then the largest box used in the 
counting is that which optimally fits within the spherical region, in order 
to avoid any spatial biasing of near-empty boxes in which there are few 
points (i.e. at the edges). This cubical subset accordingly contains a reduced 
fraction of the total library population as listed in the associated Tables, but 
still provides for a decent statistical sampling of the distribution. Figures 120 
demonstrate the box counting regression with associated confidence levels for 
the q — 0, yielding the appropriate slopes. Note that u q — > oo" cannot be 
numerically realized, so the associated values cited herein correspond to the 
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values of D q for q > 60, at which stage the estimates are observed to have 
reached a relatively stable value. 

There is a fairly narrow range in the dimensions for each library, es- 
sentially yielding a possible range of dimensions between Dq = 2.5 — 2.8, 
depending on the overall size of the sample and the choice of fit points, with 
correlation dimensions D 2 = 2.4 — 2.6 (see also Figures fTTTfTHI for the gen- 
eral form of the D q spectra). The results might suggest that the positively- 
and negatively-curved spaces yield slightly higher box dimensions than the 
flat case, although to within the cited fit error no definite determination 
may be made. In fact, to the accuracy of the fits, the Tables are virtually 
indist inguishable . 

Note that D^^ tends to approach a value near 2 in every case. This 
result has a simple geometric interpretation, and in fact provides a definite 
signature for the PSC mechanism. Since the D q values for large q can be 
interpreted as local scaling dimensions, this implies the dense regions are 
effectively two-dimensional structures. This will be further discussed in Sec- 
tion^ in comparison with reported multifractal indices for observational and 
simulated data. 

Furthermore, based on the relative consistency of the dimensions, one 
could conclude that these are signatures of the construction algorithm, and 
independent of the space in which they are built. Note that the fractal 
dimension of two-dimensional Apollonian Packing is reported to be in the 
range Dp ~ 1.31 j21], thus perhaps the D ~ 2.6 — 2.7 is a signature of 
the 3-dimensional packing. This figure should be compared with that of 
Reference [213, which reports a value of Dp ~ 2.4 for 3-D Apollonian Pack- 
ing (although the author of notes this to be "crude"). This notion of 
construction mechanism identification will be further discussed later in this 
paper. 

Based on conclusions drawn out in Reference [TO], the dimensions in 
Tables El and E| are obtained via box counting in the flat projective sub- 
space. This projection is justified, however, since the overall angular extent 
is 2co> max = Ti / 5 (see Section I5~T]) . which constitutes a rather small portion 
of the hyper/pseudospheres. It has been shown in ^U] that curvature effects 
become minimal at about this angular extent. So, in some sense, one cannot 
expect to see any signs of curvature for such a limited library. 

It should be noted that since the libraries represent only the first-level 
in the recursion, the depth ranges only to less than 3 orders of magnitude, 
beyond which point the scale approaches the size of the smallest spheres 
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and/or the mean interparticle distance. Adding a second level of recursion 
could help to better define the overall dimension, although this can lead to 
significant increases in the computation time required to perform the oper- 
ation. Future faster processors combined with efficiently-written algorithms 
will certainly be able to handle such a task. This work is currently underway 
by C. C. D. et al. 

From a formulaic approach, it is interesting to note that homogeneity 
is preserved on the first scales of iteration in the box counting procedure. 
That is, cubes of side ~ d, d/2, d/A (where d ~ n/10 is the radius of the set) 
completely cover the distribution, giving counts of 8, 64, and (approximately) 
256 respectively. This "space-filling" quality of the packings is consistent 
with the design of the Swiss Cheese cosmology (i.e. intended to satisfy the 
Cosmological Principle on the largest of scales). 

As a quick check, these statistics are in relative agreement with fractal 
dimension calculations via alternative codes or methods. A shareware piece 
of software known as fd3 was obtained |27j . which calculates the box di- 
mension of an N-dimensional array of an arbitrary number of coordinates. 
The corresponding box-counting algorithm is based on the method of refer- 
ence [2E]. The fd3 software rather consistently pegs the dimension of the 
packings at D$ ~ 2.5 — 2.7, with mild variation in the correlation dimension 
D 2 ~ 2.4 — 2.6, in rather good agreement with the above results. 

4 Comparison to Observational and Simulated 
Data 

The results of the previous section suggest that the reported D F ~ 2 of large- 
scale clustering (Table cannot be reproduced by the PSC models with 
any level of certainty. This could imply either that the PSC model under 
consideration is limited in its predictive power, or equally that the cited 
Dp ~ 2 is incorrect. It may be that sample size plays a critical role here: 
the actual surveys have considerably smaller counts than the PSC models 
considered herein. Note that smaller models showed lower fractal dimensions, 
so while this could indicate different scaling estimates at lower ranges, it does 
not support the extension to larger scales. It is not unreasonable, then, to 
expect different scaling behavior over different ranges. 

Various multifractal analyses of the available redshift catalogs have been 
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performed, but unfortunately, there is a relative lack of these clustering anal- 
yses in the literature, as compared to the more traditional correlation func- 
tion (single fractal dimension) investigations. However, the majority of those 
which can be found tend to indicate commonalities in the data (with the ex- 
ception of studies such as those discussed previously in [11J, although this 
was not a full multifractal analysis). Recall that values of D obtained from 
the standard correlation analyses should generally correspond to q = 2 (not 
q = 0), although in some cases this tends to be a point of confusion among 
some authors and the related interpretations of analysis data. Analyses which 
are purely monofractal (correlation or conditional density) quote only a sin- 
gle exponent with no regard for the g-index, which makes for sometimes 
ambiguous cross-method comparisons. It was first noted in [2H] that from a 
multifractal analysis of the CfA cluster, there existed a mismatch between 
the box dimension (Do ~ 2.1) and the widely reported correlation dimen- 
sion (D 2 ~ 1.3), indicating that the large-scale clustering was inadequately 
described by a single scaling dimension. 

Further investigations of CfAl were reported in E] in which it was 
determined that the spectrum of generalized dimensions ranged from D ~ 2 
to Doc ~ 0.6 — 0.7. More recently, multifractal dimensions were extracted 
for the CfA2 survey jHOJ [3T|, as well as the Las Campanas survey (see 
references therein for details on each catalog). Again, a similar trend in 
dimensions was observed, although instead citing correlation dimensions with 
values D 2 ~ 2 for r ~ 5 — 30 Mpc, and strong local clustering behavior of 
-Doo ~ 1 (noting that the discrepancy in D 2 values from those reported 
previously may be due to the spatial extent of the catalog considered [31]). 
It was further noted in [31] that the multifractal behavior was only observed 
over the aforementioned distance scales, with D ~ 3 beyond (i.e. a transition 
to homogeneity ensues, again contrary to "no transition" conclusions of [B]). 

Similarly, cold dark matter (CDM) N-body simulations of gravitational 
collapse have also been studied, with essentially similar conclusions being 
drawn. These are of particular interest due to their natural connection to 
hierarchical clustering growth from small initial mass/density perturbations 
in the early Universe. The main difference, however, is that while observa- 
tional data tend to show "fluid" multifractal structure, the latter simulations 
seem to display moreover a "bi-fractal" scaling [321 EBl EUl ES] • Reference [32] 
presents a comprehensive set of dimension estimates for such simulations at 
various epochs (i.e. varying cosmological scale factor), as well for increas- 
ing biasing levels (the biasing level implies that clusters will only form for 
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peak densities above a certain threshold value, generally a multiple of the 
standard deviation of the mass fluctuation jM])- The authors demonstrate 
that any number of box dimensions between 1.2 — 2.7 are allowable, de- 
pending on the epoch and biasing levels (i.e. sensitivity to higher density 
fluctuations). Higher dimensions are observed at earlier epochs (when the 
matter distribution is still predominantly homogeneous), with the extreme 
lower dimensions occurring for high biasing. For all simulations considered, 
however, a consistent signature of the evolution is D q — > 1 for larger q, due to 
the gravitational in the overdense clustering regions (i.e. non-linear density 
perturbation regime; see e.g. [21 EH! for a review). 

References 135] support this claim for behavior in various CDM N- 
body simulations, where the high dimension D ~ 2 corresponds to larger- 
scale correlations, with D ~ 1 behavior at smaller scales (and no interim 
scaling). This motivates debate as to how well even the N-body simula- 
tions reproduce observation. In reference [33], which studies the reliability 
of various fractal dimension estimation methods, the authors discuss several 
random process simulations which reproduce both the multifractal spectrum 
of observation, as well as the bi-fractal nature of such N-body simulations. 

Thus, in addition to the mismatch between low q > dimensions, a fur- 
ther striking difference between the PSC models and the above cited results 
from both observational and simulated data are in the D q values for q — > oo. 
This is attributed again to strong gravitational clustering effects at small 
scales which dominate the densest portions of the structures ]33J |32] , whereas 
the dimensions obtained in the packing library analyses herein tend to indi- 
cate that D q ~ 2 for large q. The discrepancy is undoubtedly due again to the 
formation mechanism of the structure in question. The N-body simulations 
(and similarly the real local clustering effects) are evolutionary and dynamic, 
and hence their overall structure depends on the integrated gravitational in- 
teractions between the individual particles coupled with the scale expansion 
of the Universe. Thus, the densest structures will be "stringy", with one 
preferred direction being singled out. This is consistent with the observation 
of galactic "void- filament" structures, i.e. linear aggregate structures which 
cluster on surfaces. 

The PSC algorithm, in contrast, creates the structure "on-the-fly", with 
no consideration of gravitational interactions. Again, its construction paradigm 
is purely geometric, based on volume-optimization, and the strongest clus- 
tering will occur by smaller spheres on the surfaces of the largest spheres. 
So, ~ 2 can be interpreted as the densest regions clustering on surfaces, 
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or having clustering behavior that has only two preferred directions. This 
is certainly a reasonable statement, since by the very nature of the packing 
routine, one should expect the smallest spheres to cluster on or near the 
surfaces of the largest inscribed spheres (one can obtain a sense of this from 
Figure [TJ. These dense regions, on the other hand, are reminiscent of the 
Great Wall and similar structures observed in such redshift surveys as CfA 
(see e.g. |37j). 

As a quick clarification, the use of the term "preferred" in this case is 
used in an opposite fashion to the usual discussions of homogeneity. Note 
that Dp = 3 would be described as structure having three "preferred" scaling 
directions (i.e. instead of no preferred direction - simply put, all directions 
are equally preferential). 

The obvious mismatch in D q spectra between PSC and observation/simulation 
is discouraging, but it is not to say that the PSC model is inherently flawed. 
Rather, these figures should be seen as accepted limitations to the proce- 
dure, considering the accuracy of the data and associated statistics. The 
PSC models implicitly preserve the Cosmological Principle and Weyl flat- 
ness, while N-body simulations actively break these symmetries (or at the 
very least, do not concern themselves with consistent maintenance of this 
condition throughout the evolutionary process). 

It is perfectly reasonable, however, to suggest that the two differing mod- 
els are complementary, rather than in opposition. PSC shows large-scale 
homogeneity, preserves the CP, and displays strong clustering in two dimen- 
sions surrounding voids. All that is really missing is the filamentary struc- 
ture. That is, the N-body simulations could always take place within any 
particular sphere of the packing. Recall that although the Weyl condition 
is satisfied on the large scale, in accordance with the FRW dust solutions, 
but locally (within the sphere) it can be violated (i.e. the vacuum is Ricci 
flat, but a non- vanishing Weyl is required to maintain overall curvature). 
Thus, the Weyl constraint provides a necessary boundary condition for the 
preservation of the CP. A hybrid of different clustering paradigms could help 
to yield a more realistic approximation to the observed structures and the- 
oretical constructs. Equivalently, if the N-body simulations alone are truly 
bi-fractal, then combined with a PSC-like structure, one can obtain a more 
realistic model with which to compare observational data. 

The fact that this discrepancy between dynamical and "instantaneous" 
formation paradigms can be readily signed through the multifractal analysis, 
and that the Universe more closely resembles the N-body data, can lend a 
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sense of relief to philosophers. That is, this is a sure signal that the Uni- 
verse has been evolving, or at the very least, was not created on-the-fly in 
the manner of the PSC! Thus, again, it would be interesting to study the 
multifractal scaling behavior of an N-body simulation nested within a PSC 
library, left to evolve over time. If this were to match observational data, 
the result could have profound impact for early Universe dynamics (not to 
mention the philosophical implications of a pre-determined structure in the 
beginning) . 

Along the same lines, perhaps the observational D ~ 2 of the galaxy 
clusters is actually a measure of this quantity, constrained by the relatively 
small sample space of candidate galaxies with respect to the estimated size 
of the Universe. Hence, the PSC could conceivably represent a "large-large" 
scale structure of the Universe, that is, much larger than the current obser- 
vational limits. This would then be consistent with the above characteristics 
of observation and theory. 

A brief mention is in order for another earlier hierarchical cluster model 
of related interest, that of Soneira and Peebles [38| . Within a volume of 
space, a sphere or radius R is inscribed, and within that sphere are placed 
n spheres of radius Rj p. In turn, within each of those spheres are placed n 
new spheres of radius Rj p 2 , and so forth. It is interesting to note the concep- 
tual similarities between this model and the SC packings, both of which are 
based on instantaneous positioning of spheres in a volume, and not on any 
time-evolution paradigm. The Soneira-Peebles model is not explicitly mul- 
tifractal, but rather a simple monofractal of dimension log(n)/ log(p). The 
corresponding multifractal spectrum is effectively flat, such that D q = Dp Vg 
for r e (R/p n ~ l ,R), where n is the recursion level This, however, bears 
little resemblance to the PSC spectra, nor for that matter to observational 
data. 

4.1 Comparison with Random/ Ordered Distributions 

Since the results do not suggest any method of discrimination between the 
three classes of models, it is important to address the question of how the 
multifractal structure of the packings compare to those of other random 
distributions of points. As a test, several toy random spherical distributions 
were populated using the C library function drand48(), which on each call 
generates a pseudorandom floating point number in the range [0, 1). Since 
curvature plays minimal or no role in the resulting spectra, for conciseness 
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comparisons will be made only to the flat libraries. 

One set of distributions is designed to be completely randomized, such 
that each (x, y, z) coordinate is generated randomly. A second, dubbed 
"linear-random" , randomly generates (r, 6, </>) coordinates in the usual range. 
A third represents a completely uniform spacing of points which fills the ap- 
propriate volume. These three test sets, along with a random packing library, 
are displayed in Figures Eland El Since no apparent distinctions can be made 
between the flat, open, or closed universe distributions, only the flat library 
SCF3 will be addressed in this section. 

The simple Box dimension (-Do) shows mild differences between each set, 
and in certain cases fails to accurately distinguish between them. However, 
the larger q values do indeed show large variations in the overall structures. 
In particular, note that the linear-random distribution gives a box dimension 
of D q ~ 2.8, but rapidly decreases for q > to a value of D q — > 1 for q — > oo. 
This exemplifies the ability of the multifractal spectrum to pick out anoma- 
lous structural qualities. Although the coordinates (r, 9, 4>) are populated 
randomly, the volume element of the space is not. For randomly selected 
(x,y,z), the volume element dV = dx dy dz is itself purely random, while 
randomly selected (r, 8, 0) coordinates will yield a radially- weighted volume 
element dV = r 2 sin 9 dr d9 dip. The overall clumping will reflect this r pref- 
erence, and hence will appear "pseudolinear" . Since each set fills the volume 
almost homogeneously, and as such "fools" the box counting algorithm into 
a misrepresentation of the entire structure. This LR set can be compared 
to the local behavior of the N-body simulation dimensions discussed in Sec- 
tion El both of which are dominated by purely radial {e.g. gravitational) 
clustering behaviors at small scales. In fact, it is interesting/amusing to note 
that for this set that the correlation dimension D 2 ~ 1.9, further suggesting 
that this test data could be considered a potential candidate for large-scale 
galaxy clustering (although purely in jest; this behavior appears more bi- 
fractal, due to the quick decrease to D q ~ 1, and most likely lacks additional 
structural features beyond the D q spectra). 

The -Dqq value signals the linearity of the densest clustered regions. The 
randomly populated (x, y, z) space (RD) is more reflective of a truly random 
sample, as indicated by the dimensionality. Note that ideally, the set should 
have D q = D = 3 Vg, but the variance can be understood to be software 
and data set limited (the same would be true of UD, which is just a variant 
of the RD). 
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4.2 Ribeiro's Tolman Swiss Cheese Cosmology 

As previously mentioned, the structure of a Swiss Cheese-like cosmology is 
a natural choice if one adheres to the existence of local inhomogeneity while 
adhering to the Cosmological Principle. Similar solutions have in fact ap- 
peared in the literature before, perhaps one of the most thorough being the 
model presented in the series of papers by Ribeiro EQ 122] ■ It is interest- 
ing to compare these results - which discuss a relativistic fractal cosmology 
- to those of the current paper. Ribeiro's model is very similar in structure 
to the present Swiss Cheese model considered herein, matching Tolman so- 
lutions with FRW dust solutions and integrating local density distributions 
along the past light cone to calculate the observed fractal dimensions. For a 
variety of classes of solutions, the author finds fractal dimensions which range 
between/}^ = 1.3 — 1.7 depending on the model type under consideration, 
(subject to the constraint of obeying the de Vaucouleurs' density power law). 
The interested reader is directed to the aforementioned citations for further 
reading. 

5 Curvature Considerations 

5.1 Evaluation Along Geodesies in k = ±1 Spaces 

Ideally, one must be careful when evaluating radial distances in non-flat ge- 
ometries. Although it has been argued that the curved geometries considered 
herein possess "essentially-flat" characteristics due to their small angular ex- 
tent, it is worth a quick check to see if there might exist any overt differences 
between the flat box estimates and those which do consider curvature effects. 
The spatial portion of the FRW metric can be written as 



where the Cartesian embedding coordinates for the corresponding hypersur- 
faces are [2] 



da 2 



dcu 2 + S k (cu) 2 dO 2 + sin 2 9 dej) 2 



(9) 
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Sk(ui) sin^) cos(</>) 
Sk(uj) sin(6 ) ) sin(0) 
Sk(u)) cos(6>) 



(10) 
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Here, T k (uj) = cos(u;) for k = +1 (cosh(u;) for k = —1), and the radius of 
curvature has been set as R = 1. 

Realistically, one must compute the geodesic distance between points in 
the embedded space, instead of the straightforward Euclidean distance. For 
the positive and negative spaces considered herein, Equation H3I is evaluated 
along angular geodesic distances r = 5, which is extracted from the usual 
inner product in the space between vectors (u 1 ,9 1 ,(f) 1 ) and (u 2 , 9 2 , (fi 2 ), 

cos S = cos(a>i) cos(u; 2 ) + sin(o;i) sin(a;2)[sin(6 l 1 ) sin(6 , 2 ) cos(0i — 4> 2 ) 

+ cos(0i)cos(0 2 )] , (11) 

for k = +1. The analogous expression for k = — 1 obtained by the usual 
replacement for Sk(u>) and the appropriate sign change, 

cosh 5 = cosh(u;i) cosh(c<j 2 ) — sinh(c<j 1 ) sinh(u;2)[sin(0i) sin(0 2 ) cos(0i — (f) 2 ) 
+ cos(0i)cos(0 2 )] , (12) 

The scaling dimensions may then be evaluated with "bottom-up" estima- 
tion techniques along the geodesic. The conditional average density T(r) of 
points in the set, defined as 

r(r) _ 1 f 1 dN l{ r) 3 (13) 

estimates the average change in number of points within a spherical shell of 
radius r = 5 (area A s (r)). The shell area A s (r) will depend on the value of 
k, however these effects are exceedingly small and can be ignored. This is re- 
ported to be a better estimation tool than the standard correlation function 
£(r) = (r/Yo) D ~ 3 — 1 (where Dp = D 2 ), which makes a priori assumptions 
of homogeneity with regard to the set under consideration (for an inhomo- 
geneous fractal object can lead to spurious results) 0. Of course, disagree- 
ments between dimension estimates from these two methods can also arise 
if the distribution is not a pure monofractal. If the set in question is frac- 
tal, then the number N(r) oc r D , and thus T(r) oc r D ~ 3 , which provides a 
good measure of the fractal dimension D (and is frequently used in the cited 
references). Evaluating (|13|) over iV ~ 2000 points per set or less, with an 
equivalent i? e ff < 0.10, one can obtain figures in good agreement with the 
cite box counting estimates. 



17 



The plots of Figures ISHTU1 show T(r) for sample libraries of each curvature, 
as compared with a base slope of 2.8. The fit itself is somewhat dependent 
on the choice of points, but all seem to suggest a similar slope. The smaller 
population sets show higher variability in the convergence of points to linear- 
ity, but this effect is smoothed out as the set grows in size. Note that the flat 
libraries tend to yield slightly smoother trends than the positive or negative 
sets. While this could be a curvature effect, it is unwise to make such an as- 
sertion without further investigation. The apparent non-linearity of each set 
below r ~ 0.01 is a recognized artifact of finiteness in volume-limited samples 
such as these, as the radial distances drop below the range of statistically- 
significant clustering |41j . Hence, the cited dimensions are obtained for the 
linearity which ensues to the right of this "peak" This method is useful as 
a self-check here, since the calculated dimensions tend to support those ob- 
tained via the box counting method. Since it cannot be used to measure 
additional multifractal structure, however, its uses have been exhausted. 

6 D q for q < 

The value of the parameter q need not be restricted to positive integers, nor 
in fact need it be restricted to integers. The D q are well-defined for all q G 71. 
Whereas the D q values for q > represent the scaling behavior of increas- 
ingly dense regions, those for q < correspond to the scaling of under-dense 
regions. Both the analyses of observational data and N-body simulations 
mentioned in Section |U suggest that D q — > 3 for q — > oo, indicating that in 
all cases the least populous or dense regions scale effectively as homogeneous 
distributions. 

Numerically, obtaining D q values of q < can be a very difficult quantity 
to estimate via box counting, since calculations becomes severely dependent 
on the finite size and population of the data set [33J. Such behavior is ob- 
served herein - the associated error estimates grow significantly for q > 0, 
particularly for values of q < —1.5 or so. For example, the box counting 
algorithm yields D q > 2.7 (0.3) for q = -2 for SCF3. Similarly, SCP5 
shows D q ~ 2.8 (0.2), and 2.7 (0.3) for SCN3. The other libraries considered 
demonstrate roughly similar behavior, although below these values of q, the 
fit errors quickly grow. In fact, in certain cases the estimated D q values begin 
to decrease for sufficiently small q, but with significantly increasing error (for 
all points fit). This could indicate that accurate determination of such gen- 
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eralized dimensions is extremely dependent on the choice of points returned 
from the calculation. 

There are several other means by which one can more accurately compute 
these values. A common method used in astrophysical analyses is the den- 
sity reconstruction algorithm (see e.g. jSJ E2] and related references), which 
determines the minimal radius r(p) around a point for which the probability 
is p of finding p N tot points. The corresponding partition function 

i Ntot 

w(r,p) = — - ]T nip)-^ ex p 1 -* , (14) 

iVtot i= i 

Recall that the multifractal spectrum {D q } is defined by the relation r = 
r (<z) = (q ~ 1) D q , thus the corresponding generalized dimensions may be 
obtained accordingly. This method is claimed to converge well for q < — 1, 
but is still applicable for a small range of positive q. 

The calculation is somewhat computationally intensive, since ideally the 
radial calculations must be made for every point in the set. Generally, a 
smaller (random) sampling of the set can be used to evaluate the partition 
function. Sample calculations via implementation of this algorithm have been 
implemented for a range of probability values p £ [0.01,0.1] (similar to the 
range used in [5]), for points whose inscribed radius is not overly large (else 
the calculation becomes saturated and useless), which generally includes up 
to several thousand points. In the cases of k = ±1, these are evaluated 
out along the geodesies (although as has been established, this complication 
likely isn't necessary). 

Figures ITTIITBI show sample D q plots for various libraries considered herein. 
The curves are shown to demonstrate the general trend in D q values, since 
there is sometimes a mild mismatch between the two methods (although 
within the box counting error). Thus, explicit points and error bars have 
been suppressed for the time being. The errors are roughly 0.1, and so the 
curve can be assumed to be within 1- or 2-cr of the actual values. Initial 
results indicate that the scaling dimensions of the packings for small q « 
approach D q ~ 3.0, as was suggested by the box counting method. The rate 
at which the values for individual sets in each class (flat, positive, negative) 
approach the limiting value show some mild variation, but like the cases for 
positive q, this behavior could simply be due to the size of the data set under 
consideration (e.g. SCF1 shows markedly lower values than the rest, perhaps 
due to its smaller population size). 
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The density reconstruction method is also applicable to (small) values 
of q > 0, and it is noted that this provides good agreement with the box 
dimension estimates of Section EJ In particular, the estimates for q ~ 2 tend 
to be closer in agreement to the fits of Tables 120 to within the cited error, 
thus providing yet another self-check for the estimates therein. 

In certain cases, the D q value does not rise much above that of D , how- 
ever it is uncertain whether or not this is an actual artifact of the data set, or 
rather calculation anomalies. This does not appear to be an explicit artifact 
of curvature, since it occurs for packings belonging to each class. It could be 
that these values of D q can serve as some variety of identification for varying 
initial conditions or the general distribution of sphere sizes (see That 
is, while q > is a measure of the strongest clustering, which can be inter- 
preted as an overabundance of local mass at small scales, q < is to a certain 
extent a measure of open space, which certainly could vary from library to 
library. It could also be some manifestation of curvature on the construction 
paradigm at the largest levels [TO]. Certainly, further investigation into the 
behavior is ideally warranted. 

Coupled with the box counting results from before, these figures seem 
to consistently suggest an asymptotic value approaching D q ^ 

—oo ~ 3. Such 

behavior has the simple interpretation that the regions of least clumping scale 
homogeneously, and thus no explicit structural information can be extracted 
from these values. Similar behavior is reportedly observed for the N-body 
simulations discussed previously, as well as the analysis of observational data 
from the diverse galaxy catalogs. Thus, differentiation of formation models 
seems to become even more ambiguous for this range of scaling. The strong 
clustering behavior for q ^> gives a much more intuitive and exploratory 
glimpse of the inherent structure qualities of the set. 



7 Biasing Counts by Luminosity / Mass 

As previously noted, the estimated scaling dimensions for the packing li- 
braries are high due to the density of points per allocated volume. Another 
issue is the unbiased nature of the counts. Simply put, the libraries are 
geometric configurations which satisfy the relativistic field equations. Each 
point is assigned a mass, which is a function of the inscribed radius within 
the constant density field, dependent on the curvature of the spatial man- 
ifold. However, to date this value has not been factored into the counting 
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scheme. 



The observed three-dimensional distributions are interpolated from mea- 
sured catalogs - what one sees is what one gets. However, it is not unreason- 
able to suggest that what one sees is not what is actually there! The notion of 
"missing mass" is a recurring theme in many cosmological endeavors, rang- 
ing from dark matter issues, to low luminosity objects. Furthermore, it is 
common-place in observational astronomy to simply excluded objects whose 
luminosity is below a given threshold, irrespective of whether or not it is 
visible. If it does not meet the selection criteria for the catalog in question, 
it is omitted altogether. Furthermore, the issue of luminosity biasing or seg- 
regation effects is one which is frequently raised in the debate surrounding 
the D ~ 2 over all length scales. For example, it has been suggested that 
luminosity segregation favors brighter galaxies at greater distances, which 
consequently implies a biasing toward stronger clustering, which could skew 
the actual measured dimension from homogeneity to some different scale (see 
[32*] and associated references therein). 

The mass-to-light ratio for galaxies is somewhat dependent on factors such 
as age, morphological type, and so forth. It seems reasonable to impose a 
mass cutoff in the packing libraries as a first-pass gauge of this effect. That is, 
can one impose constraints on the range of lower mass cutoffs from which one 
can extract a dimension of D ~ 2 for the packing libraries? For most galaxies 
of a particular class, this ratio is roughly a constant, i.e. M(L) oc L 13 , /3 ~ 1. 
For simplicity, such a relation will be assumed for the analysis herein, and 
thus the results can provide a good "boundary" for what one might expect in 
realistic observational circumstances. Since the initial matter field is assumed 
to be of constant density po, the mass Mj of the i th sphere may be calculated 
by the usual relations 



where Ri is the inscribed (angular) radius of the sphere, for flat ("T5*|). positive 
(HHJ), and negative curvatures (JT7J). The latter expressions are obtained by 
integrating the FRW volume element for the respective values of k. Since po 
is arbitrary, one can set p Q = 1 WLOG. 

Figure Hoi shows the mass-reduced D q spectra for various reduced libraries, 
using a mass cutoff approximately 0.01% that of the largest sphere mass. In 
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terms of the assumed mass-to-light ratio, this logically amounts to a similar 
scaling for the object's luminosity. The reduced data sets contain roughly 
10% of the original number, but still accounts for over 90% of the total 
enclosed mass. Similar results are obtained for the other packing libraries 
using a similar cutoff. This indicates that the majority of the spheres are 
relatively small in mass, and thus by the mass-to-light ratio, would most 
likely be much fainter relative to the remaining points. Hence, the smaller 
masses serve to "smooth out" the overall matter distribution, leading to 
the saturation of homogeneity observed in the upper-most levels of the box 
counts. 

Note the somewhat wide variation in limiting values of D q for negative q. 
Although the data confirms sparser clustering in these regimes, the asymp- 
totic value which D q assumes tends to be somewhat dependent on the library. 
Initial results seem to suggest that there is no preferential pattern for spe- 
cific geometries, with both possessing slow and fast approaches to a higher 
limiting dimension. It may be the discussed previously, that further 

investigation into the q < regime could potentially shed light on the as- 
sociated geometry, although a better understanding of the sensitivity of the 
method to population sizes is required. The authors of [32| similarly sug- 
gest that q < values could help discriminate between differing initial power 
spectrum perturbation conditions for N-body simulations. 

Table displays the results for the associated cutoffs for each the libraries 
of Tables I2E1 Again, note that there is no explicit signature variation in the 
estimates between curvature cases. The SCP libraries show slightly higher 
Dq and lower Dqo values than the SCF and SCN packings, albeit all equal 
to within the associated error. The same mass cutoff has been applied to 
each case, even though there is mild variation in the overall masses for each 
library (depending on the packing number). However, this is most likely a 
small consideration which does not significantly affect the end result. This 
variation could also be due to volume limitation constraints which exist in 
the closed manifold, but not in the flat and open cases. 

Thus, it would appear that the use of luminosity cutoff can potentially 
bring the box and correlation dimensions of the libraries in closer agreement 
to the reported values of reference jH] (albeit with larger fit error), which 
suggests that such a biasing mechanism could help to rationalize the afore- 
mentioned discrepancy between the PSC models and observation (if correct). 
Such luminosity biasing in Abell clusters was discussed in detail by Bahcall 
and Soneira, who showed a mismatch in the galaxy and cluster spatial corre- 
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lation lengths of at least a factor of 5 [3D] as well as a strong dependence of 
the correlation function on cluster richness. The more recent analysis of the 
SSRS2 redshift survey presented in [10] further supports the conclusion that 
there exists a strong connection between statistical clustering and luminosity, 
in particular weighted toward bright galaxies (e.g. M < —21) over fainter 
ones. The authors comment that these results are largely inconsistent with 
current theoretical models, however this paper indicates that the PSC could 
offer a resolution to this observed effect over other clustering paradigms. 

Furthermore, this can help support the notion that observed galaxies 
(luminous matter) can form a fractal distribution (subset) within a largely 
homogeneous matter distribution, much the same way the distribution of 
mountain peaks is fractal, while the Earth itself is largely spherical |%2*] . 
Observationally, similar cutoff biasing has been recorded, and that the as- 
sociated fractal dimensions tend to rise for a decreasing galaxy luminosity 
threshold |H1 E2]. Thus, the ideas presented herein are consistent with re- 
ported astrophysical procedures, and could always be used to help provide 
insight into similar "missing mass" investigations. As with the fully unbiased 
model, however, the values of for all libraries considered tend to clus- 
ter around -Dqo = 2, which again does not match observational and N-body 
simulated models (whose ~ 1, as previously discussed). 

However, some caution must be exercised in interpreting such results. 
This trend may be moreover a statistical manifestation influenced by a dif- 
ferent linear fit choice than an actual structural change in the set. In fact, 
by neglecting the smallest box size count in the fits for the higher cutoff 
biased sets, the Dq values can be shown to rise to ~ 2.4 — 2.5 with smaller 
fit error, suggesting that the dimension is potentially higher than indicated 
in Table (although still lower than the estimated D ~ 2.7 for the whole 
library). There is also a marginally larger range between the differing di- 
mension estimate techniques for the mass-biased libraries. For example, the 
density reconstruction method tend to yield values of D q ~ 2.3 — 2.5 or so 
for q ~ (see the trends of Figure IT3j). indicating potential limitations of the 
box counting technique. 

Figure Uni demonstrates the trend in dependence of the D q values on cutoff 
size, both for fixed number of iterations for all cutoff levels, as well as adjusted 
numbers for smaller set populations. The estimates remain largely unaffected 
until a cutoff of > 0.001% is reached (leaving about 20% of the spheres and 
96% of the total contained mass), after which there is a relatively quick 
drop in D with respect to for the fixed number of iterations. Whereas, 
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adjusting the number of box levels (reducing by one) helps to correct the 
curve to a higher estimate. Beyond about 0.04%, the dimensions approach 
roughly the same value D < 2.0, suggesting a type of monofractal behavior. 
At this point the population of the set has been reduced to 2% of the original, 
thus the statistics could start to become skewed by fmiteness effects due to 
a larger average point separation. 

The box counting technique has at times been criticized in the literature 
for being too sensitive to discreteness effects, especially when the population 
size is small (causing spurious results or underestimates of the actual scal- 
ing dimensions) [3U 1123 • In other cases, box counting is hailed as a quite 
robust and stable method for estimation of the associated statistics, and 
furthermore the density reconstruction method is claimed to give dimension 
overestimates for limited size sets [33] ■ Thus, there is a certain level of dis- 
agreement and confusion in the community with respect to the utility of any 
of these methods. 

So, when few data points are available such as in this situation, it is 
probably best to take some variety of average dimension as calculated by 
differing methods. Although the measured D Q dimensions do drop when 
lower masses are discounted, the actual magnitude of the drop seems sensitive 
to the size of the data set (and measurement method). Whether or not 
this represents a skew in the measured dimension as a result of luminosity 
biasing is somewhat unclear. Thus, such biasing results are more susceptible 
to statistical anomalies (highly dependent on the choice of points in the box 
count fit) than are the full models, and care must be taken in their use. It 
is, however, worth noting that the values remain relatively consistent 
at ~ 2, indicating a certain robustness of the data. This is a definite 
signature of the packing structure. 

From a physical point of view, one must recall that the PSC models 
are merely weighted points. As mentioned in a previous section, one could 
imagine replacing each point with another type of set, such as an evolutionary 
N-body simulation, which could then conceivably yield finer clustering effects 
(and potentially yield a ~ 1 behavior). If the dimensions of Table are 
actually statistical properties of the luminosity-biased sets, then it could 
constitute a possible solution to reconciling the mismatch with observational 
data. 
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8 Future Considerations 



The basic premise of the packing algorithm is to ensure preservation of the 
Cosmological Principle via maintenance of Weyl flatness. Although this is 
done by compressing the matter contained within a spherical region, there 
is a priori no reason for this choice (apart from the motivation of gravita- 
tional collapse). Suppose instead the matter were expanded to lie along the 
spherical shell described by the inscription. This would create a thin sheet 
of matter, consistent with the two-dimensional structures complemented by 
voids, evident both from observation as well as the multifractal analysis. The 
matter could then be allowed to coalesce by some gravitational mechanism 
on the sheets, yielding the appropriate linear structures. Alternatively, the 
packing algorithm could randomly choose between collapse and expansion 
for the spheres, and in the case of a collapse, the interior matter is allowed 
to cluster via N-body behavior. Whether or not such models are physically 
realizable is unknown, but they could nevertheless provide yet another model 
from which to study and reconcile the observed large-scale structures in the 
Universe, and as mentioned, stand to uphold the Cosmological Principle on 
many levels. 

Furthermore, the use of the monofractal dimension can frequently short- 
change the characterization of a clustering set, since it is possible to have 
vastly differing structures which possess the same base dimension Dp = D . 
Although the multifractal spectrum can help differentiate such situations, 
a different consideration is the lacunarity of the set. This measure, often 
associated with the texture of a fractal, provides an estimate of the "voidness" 
(rather than "dumpiness") (see e.g. [TJ |2Hl EH EH] ) • Such an investigation is 
currently underway by J. R. M.. 
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survey 


D F 


Approx. Size 


LiAl 


1.7 (0.2) 


1800 


CfA2 


~ 2 


11000 


SSRS1 


2.0 (0.1) 


1700 


SSRS2 


~ 2 


3600 


LEDA 


2.1 (0.2) 


75000 


IRAS 1.2/2 Jy 


2.2 (0.2) 


5000 


Perseus-Pisces 


~ 2.1 


3300 


ESP 


1.8 (0.2) 


3600 


Las Campanas (LCRS) 


2.2 (0.2) 


25000 



Table 1: Galaxy fractal dimension calculations for various redshift surveys 
(compiled from jHj). 
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Library (pts) D D 2 D 



SCF1 (22076) 
SCF2 (42325) 
SCF3 (49797) 
SCF4 (73936) 
SCF5 (82788) 



2.5 (0.1) 

2.6 (0.1) 

2.7 (0.1) 

2.7 (0.1) 

2.8 (0.1) 



2.4 (0.1) 
2.6 (0.1) 
2.6 (0.1) 
2.6 (0.1) 
2.6 (0.1) 



2.2 (0.1) 

2.2 (0.1) 

2.2 (< 0.1) 

2.3 (0.1) 

2.3 (< 0.1) 



Table 2: Calculated dimensions for flat space packings. Cube size ranges 
from homogeneous saturation at largest scales (8 cubes) to roughly 2-4 points 
per cube (ranging from scales (~ n/10 — 7r/300)). 
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Library (pts) D D 2 D, 



SCP1 (31904) 
SCP2 (35268) 
SCP3 (47140) 
SCP4 (54966) 
SCP5 (80437) 



2.6 (0.1) 

2.7 (0.1) 
2.7 (0.1) 

2.7 (0.1) 

2.8 (0.1) 



2.5 (0.1) 

2.6 (0.1) 
2.6 (0.1) 
2.6 (0.1) 
2.6 (0.1) 



2.2 (0.1) 
2.2 (0.1) 
2.2 (0.1) 

2.2 (0.1) 

2.3 (0.1) 



Table 3: Calculated dimensions for positively-curved space packings. 
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Library D D 2 D, 



SCN1 (32965) 
SCN2 (49343) 
SCN3 (60245) 
SCN4 (65899) 
SCN5 (83863) 



2.6 (0.1) 

2.7 (0.1) 
2.7 (0.1) 

2.7 (0.1) 

2.8 (0.1) 



2.5 (0.1) 

2.6 (0.1) 
2.6 (0.1) 
2.6 (0.1) 
2.6 (0.1) 



2.1 (0.1) 

2.2 (0.1) 

2.3 (0.1) 
2.3 (0.1) 
2.3 (0.1) 



Table 4: Calculated dimensions for negatively-curved space packings. 



32 



Lib. 


# (% tot) 


7o mass 


D 


D 2 


D 


oo 


SCF1 


4526 


(20.5) 


96.0 


2.1 


(0.3) 


2.1 


(0.2) 


1.9 


[0.1) 


bLr 2 


5074 


(12.0) 


93.6 


2.1 


(0.3) 


2.1 


(0.2) 


1.9 


v 0.1) 


SCF3 


5218 


(10.5) 


92.8 


2.2 


(0.3) 


2.2 


(0.2) 


1.9 


[0.1) 


SCF4 


5506 


(7.4) 


91.1 


2.2 


(0.3) 


2.2 


(0.2) 


1.9 


[0.1) 


SCF5 


5599 


(6.7) 


90.6 


2.2 


(0.3) 


2.2 


(0.2) 


1.9 


[0.1) 


SCP1 


4744 


(14.9) 


94.6 


2.1 


(0.3) 


2.1 


(0.2) 


1.8 


[0.2) 




5848 


(16.6) 


no o 
93.3 


2.3 


(0.3) 


2.2 


(0.2) 


o n 
Z.U 


v 0.1 ) 


SCP3 


5108 


(10.8) 


92.9 


2.2 


(0.3) 


2.2 


(0.2) 


2.0 


[0.1) 


SCP4 


5233 


(9.5) 


92.4 


2.2 


(0.3) 


2.2 


(0.2) 


1.9 


[0.1) 


SCP5 


5542 


(6.9) 


90.6 


2.2 


(0.3) 


2.2 


(0.2) 


1.8 


[0.2) 


SCN1 


4993 


(15.1) 


94.7 


2.1 


(0.3) 


2.1 


(0.2) 


1.9 


[0.1) 


SCN2 


5294 


(10.7) 


93.0 


2.2 


(0.3) 


2.2 


(0.2) 


1.9 


[0.1) 


SCN3 


5409 


(9.0) 


92.1 


2.2 


(0.3) 


2.2 


(0.2) 


1.9 


[0.1) 


SCN4 


5536 


(8.4) 


91.7 


2.2 


(0.3) 


2.2 


(0.2) 


1.9 


[0.1) 


SCN5 


5837 


(7.0) 


90.6 


2.2 


(0.3) 


2.2 


(0.2) 


2.0 


[0.2) 



Table 5: Mass-reduced box counts for all libraries, showing box, correla- 
tion, and limiting Dg^^ values. Lower D values with higher errors may be 
artificially lower due to finiteness effects. 
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Figure 1: Sample projected two-dimensional visualization of packing library, 
containing approximately 35000 spheres. Surfaces represent the inscribed 
radius, and not the actual contained matter. 
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Figure 2: Box counting (£>„) fit for SCF3, yielding D = 2.7 (0.1). 
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Figure 3: Box counting (D ) fit for SCP3 yielding D = 2.7 (0.1). 
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Figure 4: Box counting (D ) fit for SCN3 yielding D = 2.7 (0.1). 
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Figure 5: Three-dimensional point distributions, including: packing library 
(SCF3), randomly-populated (x,y,z) coordinates (RD). Each set possesses 
between 50000-60000 points. 
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Figure 6: Three-dimensional point distributions, including linear-random 
(LR; top), and uniform distribution (UD; bottom). Each set possesses be- 
tween 50000-60000 points. The LR set is comprised of randomly-populated 
(r, 9, 0) points, in contrast to the random population of (x, y, z) coordinates. 
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Figure 7: D q spectra (q > 0) for distributions of Figure 15161 All sets show 
high D values, with the RD and UD sets (top two curves) maintaining 
D ~ 3. Explicit structural differences can be seen in the SC and LR sets 
(bottom two curves) for large q, signifying strong inhomogeneous clustering 
behaviors. 
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h 6 



3 I 

0.001 



SCF1 
SCF2 
SCF3 
SCF4 
SCF5 




++++++ + 



Figure 8: Conditional density T(r) for flat libraries. Sample fit (dashed line) 
has slope D — 3 = —0.2 , which implies D = 2.8. Note that each library has 
been shifted vertically to eliminate overlap for easier viewing. 
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Figure 9: Conditional density T(r) for positive libraries. Sample fit (dashed 
line) has slope D — 3 = —0.2 , which implies D = 2.8. 
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Figure 10: Conditional density T(r) for negative libraries. Sample fit (dashed 
line) has slope D — 3 = —0.2 , which implies D = 2.8. 




Figure 11: Bezier curve fit of full D q spectra for flat packing libraries, show- 
ing D g 3 for small q. q < values calculated by density reconstruction; 
q > by box counting. 
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Figure 12: Bezier curve fit of full D q spectra for positive packing libraries. 



o° 2.5 




Figure 13: Bezier curve fit of full D q spectra for negative curvature packing 
libraries. 
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Figure 14: Library SCF3 with effective mass cutoff, containing 10% of 
spheres representing approximately 90% of total contained mass. 
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Figure 15: Bezier curve fit of full D q spectra for mass- reduced libraries, show- 
ing D q 3 for small q. q < values calculated by density reconstruction; 
q > by box counting. 
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1.6 

1e-05 0.0001 0.001 
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Figure 16: Variation of Dq (A, C) and Doo (B, D) as a function of cutoff 
(% of largest mass) for fixed number of box level iterations (A, B) and "last- 
point-removed" (C, D). True fractal dimension is bounded by the two curves. 
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